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Assessing the role of interference in natural and artificial quantum dyanamical processes is a crucial task 
in quantum information theory. To this aim, an appopriate formalism is provided by the decoherent histories 
framework. While this approach has been deeply explored from different theoretical perspectives, it still lacks 
of a comprehensive set of tools able to concisely quantify the amount of coherence developed by a given dy¬ 
namics. In this paper we introduce and test different measures of the (average) coherence present in dissipative 
(Markovian) quantum evolutions, at various time scales and for different levels of environmentally induced de¬ 
coherence. In order to show the effectiveness of the introduced tools, we apply them to a paradigmatic quantum 
process where the role of coherence is being hotly debated: exciton transport in photosynthetic complexes. 
To spot out the essential features that may determine the performance of the transport we focus on a relevant 
trimeric subunit of the FMO complex and we use a simplified (Haken-Strobl) model for the system-bath interac¬ 
tion. Our analysis illustrates how the high efficiency of environmentally assisted transport can be traced back to 
a quantum recoil avoiding effect on the exciton dynamics, that preserves and sustains the benefits of the initial 
fast quantum delocalization of the exciton over the network. Indeed, for intermediate levels of decoherence, the 
bath is seen to selectively kill the negative interference between different exciton pathways, while retaining the 
initial positive one. The concepts and tools here developed show how the decoherent histories approach can be 
used to quantify the relation between coherence and efficiency in quantum dynamical processes. 
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I. INTRODUCTION 


Coherence is ultimately the most distinctive feature of 
quantum systems. Finding a proper measure of the coherence 
present at different times scales in a quantum dynamical sys¬ 
tem is the first essential step for assessing the role of quan¬ 
tum interference in natural and artificial processes. This is of 
particular relevance for those quantum evolutions in which in¬ 
formation or energy are transformed and transferred in order 
to achieve a given task with high efficiency. In this context 
the relevant questions are: how much and what kind of co¬ 
herence is created vs destroyed by the dynamical evolution? 
How does coherence determine/enhance the performance of 
the given process? These are in general difficult questions and 
to be answered they require an appropriate and sufficiently 
comprehensive framework. A general and fundamental for¬ 
malism to describe quantum interference is provided by the 
decoherent histories (DH) approach to quantum mechanics. 
DH have mainly found applications to foundational issues of 
quantum mechanics such as the formulation of a consistent 
framework to describe closed quantum systems, the emer¬ 
gence of classical mechanics from a quantum substrate, the 
solution of quantum paradoxes, decoherence theory, quantum 
probabilities M- However, DH can also be a systematic 
tool for quantifying interference in quantum processes, and 
discussing its relevance therein. Indeed, DH provide a precise 
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mathematical formalization of interference by means of the 
the so called decoherence matrix V. The latter is built on the 
elementary notion of histories and allows one to describe the 
quantum features vs the classical ones in terns of interference 
between histories, or pathways if one resorts to the mental 
picture of the double slit experiments. It is however difficult 
to quantify in a compact and meaningful way the content of V 
and its implications for the dynamics of specific systems. Our 
first main goal is therefore to define and test appropriate mea¬ 
sures allowing for the investigation of how interference can 
determine the performance of a given quantum information 
processing task. Starting from V and by its sub-blocks we de¬ 
fine different functionals. In particular, we introduce a global 
measure of coherence C able to describe the coherence con¬ 
tent of a general quantum evolution at its various time scales; 
an average (over different time-scales) measure of coherence; 
and average mesure of interference between histories leading 
to a specific output. 

While the tools we introduce are of general interest and ap¬ 
plication, in order to test them we apply them to a specific but 
relevant instance of quantum dynamics taking place in photo¬ 
synthetic membranes of bacteria and plants: quantum energy 
transport. Here the basic common mechanism is the follow¬ 
ing: a quantum excitation is first captured by the system and 
then migrates through a network of sites (chromophores) to¬ 
wards a target site, e.g., a reaction center, where the energy 
is transformed and used to trigger further chemical reactions. 
There is now an emerging consensus that efficient transport 
in natural and biologically-inspired artificial light-harvesting 
systems builds on a finely tuned balance of quantum coher¬ 
ence and decoherence caused by environmental noise ll9l- [l3|] . 
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a phenomenon known as environment-assisted quantum trans¬ 
port (ENAQT). This paradigm has emerged with clarity in re¬ 
cent years, as modern spectroscopic techniques first suggested 
that exciton transport within photosynthetic complexes might 
be coherent over appreciable timescales S. Indeed, a grow¬ 
ing number of experiments has provided solid evidence that 
coherent dynamics occurs even at room temperature for un¬ 
usually long timescales (of the order of 100/s) Sd. Ef¬ 
forts to describe these systems have led to general models of 
ENAQT EliMIll, depicting the complex interplay of three 
key factors: coherent motion, i.e., quantum delocalization of 
the excitation over different sites, environmental decoherence, 
and localization caused by a disordered energy landscape. So 
far, the presence of coherence in light-harvesting systems has 
been qualitatively associated to the observation of distinc¬ 
tive ‘quantum features’. Originally, coherence was identified 
with ‘quantum wavelike’ behavior as reflected by quantum 
beats in the dynamics of chromophore populations within a 
photosynthetic complex. Later works, employing quantum- 
information concepts and techniques, have switched attention 
towards quantum correlations between chromophores, in par¬ 
ticular quantum entam ^ementil^-Id- Besides being open to 
criticism (see, e.g., 11^ l3l|] ). these approaches do not pro¬ 
vide direct quantitative measures of coherence in the presence 
of noise. Therefore, in what follows, we shall apply the novel 
tools based on DH to a simple yet fundamental model of quan¬ 
tum energy transfer. We will focus on a relevant trimeric sub¬ 
unit of the Eenna-Matthews-Olson (EMO) complex, the first 
pigment-protein complex to be structurally characterized Q. 
The trimer is virtually the simplest paradigmatic model retain¬ 
ing the basic charcteristics of a disordered transfer network 
and it can also be conceived as an essential building block of 
larger networks. Eor simplicity, we will use the well-known 
Haken-Strobl model lE^ to describe the interplay between 
Hamiltonian and dephasing dynamics. While the model is 
an oversimplified description of the actual dynamics taking 
place in real systems, it allows to spot out the essential fea¬ 
tures that may determine the high efficiency of the transport. 
We shall initially focus on a new coherence measure C, based 
on the decoherence matrix, and characterize its behavior veri¬ 
fying that it can consistently identify the bases and timescales 
over which quantum coherent phenomena are present during 
the evolution of the system. We shall then show how the av¬ 
erage coherence exhibited on those time scales can be con¬ 
nected with the delocalization process. A more detailed anal¬ 
ysis will be aimed at distinguishing between constructive and 
destructive interference affecting the histories ending at the 
site where the excitation exits the photosynthetic structure. 
By using the decoherence functional, we will show that the 
beneficial role of dephasing for the transport efficiency lies in 
a selective suppression of destructive interference, a fact that 
has been systematically suggested in the literature, but never 
expressed within a general and comprehensive framework that 
allows the quantitative evaluation of coherence and its effects. 

The application of the introduced tools and methods based 
on DH to a simple yet paradigmatic system shows how one 
can properly quantify the coherence content of a complex 
quantum dynamics and elucidate the role of coherence in de¬ 


termining the overall efficiency of the process 1^ . 

The paper is structured as follows. In Section|n]we review 
the basic decoherent histories formalism. In Section|ni]we de¬ 
fine the measure of coherence C and describe its meaning and 
properties. In Section|IV]we first introduce the used model for 
describing the energy transport in the selected trimeric com¬ 
plex. We then discuss the coherence properties of the exci- 
tonic transport: By means of the appropriate measures based 
on the decoherent histories formalism we identify the essential 
features that may determine the high efficiency of the trans¬ 
port. In Section|V]we briefly discuss how to extend our results 
to the whole EMO complex. In Section|VT]we summarize our 
results an draw our conclusions. 


II. DECOHERENT HI S TORIES 

The formalism of decoherent (or consistent) histories was 
developed in slightly different flavors by Griffiths 11112] > Gell- 
Mann SHtl, Hartlei^ and Omnes ll]. DH provide a con¬ 
sistent formulation of quantum mechanics where probabili¬ 
ties of measurement outcomes are replaced by probabilities of 
histones. In this formulation, external measurement appara¬ 
tuses are not needed, and then one does not need to postulate 
a “classical domain" of observers. As a consequence, quan¬ 
tum mechanics becomes a theory that allows the calculation 
of probabilities of sequences of events within any closed sys¬ 
tem, including the whole universe, without the necessity of 
invoking postulates about the role of measurement. In this 
framework, the “classical domain" can be seen to emerge as 
the description of the system becomes more and more coarse¬ 
grained. 

The idea of ‘histories’ stems from Eeynman’s ‘sum-over- 
histories’ formulation of quantum mechanics. As is known, 
any amplitude {iljf\U{tf — fo)|'*/’i) between an initial and a 
final state can be expressed as a sum over paths, or histories: 
Upon inserting the identity decomposition I = \j){j\ = 

Pj differerent times fi... fjv we get 

- io)|V'i) =(V’/|C^(U -tN)Y^PjNU{tN -tN-l) ■.. 

3N 

•. • U(f2 - ii) ^ - to)\^i) = 

31 

~ "y - (V’/I-Pjw (f^r) . . . Eji (tl)|'/i) 

31...3N 

where we use the Heisenberg notation Pj(t) = U\t — 
to)PjU{t — to). Thus the total amplitude ('(/;/|(7(f/ — fo)|/’i) 
is decomposed as a sum of amplitudes, each one correspond¬ 
ing to a different history identified by a sequence of projectors 
P P 

^JN ■ ■ ■ 

The decoherent histories formalism assumes that histories 
are the fundamental objects of quantum theory and gives a 
prescription to attribute probabilities to (sets of) histories. 
A history is defined as a sequence of projectors at times 
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ti < ■ ■ ■ < tN- Probabilites can be assigned within exhaus¬ 
tive sets of exclusive histories, i.e., sets of histories Sjy — 
{f 1 ,..., f AT, Pji,..., } where subscripts ji,..., Jn label 

different alternatives at times ti,... ,tN- Histories are ex¬ 
haustive and exclusive in the sense that the projectors at each 
time satisfy relations of orthogonality PjPk = SjkPj, and 
completeness, Pj = I. In other words, the projectors 
Pj define a projective measurement. Within a specified set, 
any history can be identified with the sequence of alternatives 
j = ji,..., jAT realized at times ti ,..., 

Different alternative histories can be grouped together with 
a procedure called coarse-graining. Starting from histories j 
and k we can define a new, coarse-grained history m = j V k 
by summing projectors for all times tf, such that je and kg 
differ; 

Pmi — if 7^ k^ 

Pmi_ 11^ Jr 

for all £ = 1,..., By iterating this procedure, one can ob¬ 
tain more and more coarse-grained histories. A special type 
of coarse-graining is the temporal coarse-graining: we group 
together histories j, k,..., 1 such that such that at some time 

we have -\- Pk^ Pi^. =1. Then the coarse-grained 
history m=jVkV--Vl contains only one-projector (equal 
to the identity) at time ti, that can be neglected and hence re¬ 
moved from the string of projectors defining the history. On 
the other hand temporal fine-graining can be implemented for 
example by allowing different alternatives at a times tk ^ 
{ti, In particular, one can create new sets of histo¬ 


ries Sn+ 1 = tN+i,Pji ,..., Pjjv, } from a 

given one Sn by adding different alternatives at time fAr+i > 
In', the sets Sn-^i are fine grained versions of the sets Sn- 
Once we specify the initial state g and the (unitary) time evo¬ 
lution U (f), we can assign any history j a weight 

W'i = Tr[Cj gCj ], with Cj = {In) ... Pp (ti) 

where we use the Heisenberg notation Pj^iti) = 
U {ti)PpU {ti)^. When the initial state is pure, 
Q = and the final projectors are one dimen¬ 
sional, Pj^ = l> this formula takes the simple 

form of a squared amplitude 

Wj = \{'fij^\Pj^_^{tN-i) ■ ■ ■ Pn{^i)\'fii)?■ ( 1 ) 

Weights cannot be interpreted as true probabilities, in general. 
Indeed, due to quantum interference between histories, the Wj 
do not behave as classical probabilities. Indeed, consider two 
exclusive histories j,k S 5 and the relative coarse-grained 
history m = j V k by: P^i = Pji + Pk^,^^. If the Wj were 
real probabilities, we would expect Wm = tt'j -I- Wk- Instead, 
what we find is 

Wm = Wj + Wk + 2Pe(Tr[CjpC'k]). 

Due to the non-classical term i?e(Tr[CjpC'^]), representing 
quantum interference between the histories j and k, the clas¬ 
sical probability-sum-rule is violated. The matrix 


Pjk = Tr[CjpC'^] = Tr[Pj„t/(fAr — In-i) ■ • ■ Pj.^U{ti)QU{ti)^Pk^ ... U{tN — (2) 


is called decoherence functional or decoherence matrix. The 
decoherence matrix can be thought of as a “density matrix 
over histories”: Its diagonal elements are the weights of his¬ 
tories and its off-diagonal elements are interferences between 
pairs of histories. The decoherence matrix has the follow¬ 
ing properties; i) it is Hermitian ii) it is semipositive definite 
iii) it is trace one iv) it is block-diagonal in the last index, 
Pjk = <5jjvfeN^jk- Weights of coarse-grained histories can be 
obtained by summing matrix entries in an n x n block of the 
decoherence matrix corresponding to the original fine-grained 
histories. For istance, the weight of history m = j V k is ob¬ 
tained by summing entries of a 2 x 2 block of the decoherence 
matrix; 

Wm = Djj -b Dkk + 22jk -b Dkj • 

A necessary and sufficient condition to guarantee that the 
probability sum rule wjvk = wj -b Wk apply within a set of 
histories is 


i?e[Djk] = 0,Vj ^ k. 

This condition is termed as weak decoherence^', the neces¬ 
sary and sufficient condition that is typically satisfied /jH/ and 
that we will adopt in the following is the stronger one termed 
as medium decoherence 

Djk = 0,Vj^k. (3) 

Medium decoherence implies weak decoherence. Any 
exhaustive and set of exclusive histories satisfying medium 
decoherence is called a decoherent set. The fundamental rule 
of DH approach is that probabilities can be assigned within 
a decoherent set, each history being assigned a probabability 
equal to its weight. If medium decoherence holds, the 
diagonal elements of the decoherence matrix can be identified 
as real probabilities for histories and we can write Djj = pj. 
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Due to property iv), if we perform a temporal coarse- 
graining over all times except the last, we obtain ‘histories’ 
with only one projection, Pj^ at the final time Im- These 
histories automatically satisfy medium decoherence: 

z i: = 

~ (fAf)P-Pjw (f at)] = ^JNkwPjN 

where= Tr[Pj^ {tN)QPjN (^w)] '^he probability that the 
system is in j m at time t m ■ Due to interference, the probability 
of being in jjv at time is not simply the sum of probabili¬ 
ties of all alternative paths leading to jn, i.e, of all alternative 
histories with final projection . In formulas. 

Pm 7 ^ 

The probability and the global interference of histories (''') 
ending in jV can be thus expressed as 

P3Nij)= X! (4) 

h,-m-i 

with T = NlS.t. Destructive interference will happen when 
'P-jN < 0, constructive interference when > 0. 

The decoherent histories formalism is consistent with 
and encompasses the model of environmentally induced 


decoherence||3l- Given a factorization of the Hilbert space into 
a subsystem of interest and the rest (environment), PL — Pis ® 
Pis, the events of a history take the form Pjij '= Pji ® t 
where Pj^ and Hj' are projectors onto Hilbert subspaces of 
Pis and Pis respectively. Histories for S alone can be ob¬ 
tained upon considering appropriate coarse-grainings over the 
degrees of freedom of the environment, such that the events 
are Pj^ (g> Ig where Is is the identity over Pis- Upon intro¬ 
ducing the time-evolution propagator ICt to as g{t) — U{t — 
to)gito)U{t - fo)''' = JCtto [£'(fo)]j we can rewrite the deco¬ 
herence matrix as: 

T^jk =Tr[Pj^^tiV tN-l[PjN-l^tN-l tiV-2[ (5) 

• ■ • to [t?o] • ■ • ]-Pfeiv-i] A:jv]' 

If the initial state is factorized g{to) = gs{to) ® gs{to), then 
the reduced density matrix p 5 (f) = Trs[g{t)] evolves accord¬ 
ing to gs{t) = ICt to 8s{to) where jC is the (non-unitary) re¬ 
duced propagator defined by 

T'i£[U{t-to)g~s{to)^gs{tQ)U\t-to)] = ^tto[8s{to)]- ( 6 ) 

If the evolution of the system and environment is Markovian, 
we can write iCt t’ = ^t-t'- As proved by Zurek ||3l, under 
the assumption of Markovianity we can rewrite the decoher¬ 
ence matrix in terms of reduced quantities alone, i.e., quanti¬ 
ties pertaining to the system only: 


J 


Pjk —TTtslPjN^tN tN-l[PjN-l^tM-l tN-oi- ■ ■ toiso] ■ ■ - jPkN-APkN]- 


That is, the model of environmentally induced decoherence 
can be obtained by applying the decoherent histories for¬ 
malism to system and environment together, and by coarse- 
graining over the degrees of freedom of the environment. 


III. THE COHERENCE MEASURE C 

The DH approach provides the most fundamental frame¬ 
work in which the transition from the quantum to the classi¬ 
cal realm can be expressed. Indeed, it is based on the most 
basic feature characterizing the quantum world: interference 
and the resulting coherence of the dynamical evolution. De¬ 
spite being a well developed field of study, the DH history ap¬ 
proach lacks for a proper global measure of the coherence pro¬ 
duced by the dynamics at the different time scales. We there¬ 
fore introduce a measure that quantifies the global amount of 
coherence within a set of histories. Assume projectors for 
all times £ = 1,... ,N are taken in a fixed basis \ej), 
Pj = \ej){ej\. Assume further that histories are composed 


by taking equally spaced times between consecutive projec¬ 
tions i.e., ti — At,... ,t IS = N At. (in other words, histories 
correspond to projections applied in the same basis and re¬ 
peated at regular times). For such a set of histories, consider 
the decoherence matrix 

where _ p-^{^]S[At).. .Pj,,{At). Take the von 

Neumann entropy of the decoherence matrix, 

h{P, N, At) = -Tr[I>(^’log AAi)] ^ (g) 

Due to coherence between histories, hjs differs from the 
‘classical-like’ Shannon entropy of history weights 

{P, N, At) = -J2 log (9) 

j 

where are the diago¬ 

nal elements of D(Af,P,At)^ j g ^ (jjg weights. The difference 
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between the two quantities is wider if off-diagonal elements 
of the decoherence matrix are bigger, i.e., if the set of histo¬ 
ries is more coherent. Let us define: 


CiP,N,At) 


At) - h{P,N, At) 
h(^){P,N, At) 


( 10 ) 


We argue that C{P^ N, At) is suitable to be used as a general 
measure of coherence within the set of histories defined by 
P, N, At. Indeed, we can readily prove the following proper¬ 
ties: 

i) 0 < C{P, N, At) < 1. C{P, N, At) < 1 is obvious. To 

prove C{P, N, At) > 0, let us define a matrix _ 

where off-diagonal entries are set to zero. Since 

j 

we obtain that the numerator of (fTOl) can be expressed as a 
quantum relative entropy: 

(P, N, At) - h{P, N, At) = 

= - log'D(^’-^’^O)] 

where h{A\\B) > 0 is the relative entropy between A and B. 

ii) C(P,7V,Af) = 0 iff i-e., 

C{P,N,At) vanishes if medium decoherence holds for the 
set of histories, since the two quantities h{N, P, At) and 

{N, P, At) coincide in this case. 

Thus C{P, N, At) is in essence a (statistical) distance be¬ 
tween the decoherence matrix V and the corresponding di¬ 
agonal matrix renormalized so that its value lies 

between 0 and 1. The greater are the off-diagonal elements 
of T>(PAf,At)^ j-jjg greater the distance. The meaning of 
C{P, N, At) can be easily understood if we use the linear en¬ 
tropy, a lower bound to the logarithmic version: 

hL{P,N,At) = 1 -Tr[(pl^’'f'’^*l)2], 

1 - h^^\p,N,At) = 


The measure introduced is well grounded on physical con¬ 
siderations. In the following we will apply it to a simple sys¬ 
tem in order to check its consistency, and later use it to char¬ 
acterize the coherence properties of the evolutions induced by 
various regimes of interaction with the environment. First, 
one has to check whether the measure properly takes into ac¬ 
count the action of the bath. In particular, if the bath is char¬ 
acterized by a decoherence time 7 “^, it is known (|| 3 l) that on 
time scales Af > 7 “^ the decoherence matrix becomes diag¬ 
onal: The probability of a history at time tj^^i can be fully 
determined by its probability at time tj^, since no interference 
can occur between different histories. Indeed, the action of the 
bath is to create a decoherent set of histories that are defined 
by a proper projection basis: the pointer basis (||3l)- There¬ 
fore, the fine-graining procedure obtained by constructing a 
set of histories Sn+i via the addition of a new complete set 
of projections in the same basis at time fjv+i = {N + l)At 
to the set Sn, should leave the coherence functional C in¬ 
variant, i.e., C{P,N + l,Af) Ri C{P,N,At). If instead 
At < 7 “^ the same fine-graining procedure should lead to 
C{P, N + 1, At) > C{P, N, At). 

Before passing to analyze a specific system, we want to fo¬ 
cus on the complexity of the evaluation of V and C. The di¬ 
mension of the decoherent matrix grows with the dimension 
d of the basis P and the number N of time instants that de¬ 
fine each history as . This exponential growth in princi¬ 
ple limits the application of the DH approach to small sys¬ 
tems. However, as for the system considered in this paper the 
computational effort is contained due to the small number of 
subsystems (chromophores) and the small dimension of the 
Hilbert space which is limited to the single-exciton manifold. 
As we shall see, by limiting the choice of A to a reasonable 
number, the analysis can be fruitfully carried even on a laptop. 


IV. TRIMER 


We now start to analyze decoherent histories in simple 
models of energy transfer comprising a small number d of 
chromophores (sites). Neglecting higher excitations, each site 
i can be in its ground |0)i or excited 11)^ state. We work in the 
single-excitation manifold, and define the site basis as 

|i) = | 0 )i ... |l)i... \0)d i = l...d, 


In this case, we obtain a ‘linear entropy’ proxy of C{N, P, At) 
as: 




hf{P,N,At) 




i.e., state |i) represents the exciton localized at site i. On-site 
energies and couplings are represented by a Hamiltonian H 
that is responsible for the unitary part of the dynamics. Inter¬ 
action with the environment is implemented by the Haken- 
Strobl model, that has been extensively used in models of 
ENAQT ifisl 1^ 1^ 1^ . The effect of the environment is 
represented by a Markovian dephasing in the site basis, ex¬ 
pressed by Lindblad terms L in the evolution, as follows: 


which is a simplified version that, by avoiding the diagonaliza- 
tion of helps containing the numerical complexity. 


g = [H, g] + ^7,[2L,pLJ - (11) 

i 
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where Li = |i)(i| are projectors onto the site basis, and 7 ^ are 
the (local) dephasing rates. Furthermore, site d can be inco¬ 
herently coupled to an exciton sink, represented by a Linblad 
term 


ktrap[‘^LtrapQL^j.g^p L^^g^pLtrapQ Q^trapk^trap] 


where Ltrap = \sink){e\ and ktrap is the trapping rate. Con¬ 
trary to other works, we neglect exciton recombination, as it 
acts on much longer timescales (~ 1 ns) than dephasing and 
trapping. 

The global evolution is Markovian and can be represented by 
means of the Liouville equation 


g = C{q) = CH{g) + + Ltrapig) ( 12 ) 

that can be simply solved by exponentiation. 


energy transfer begins, while the third chromophore is the site 
from which the excitation leaves the complex. The Hamilto¬ 
nian of the trimeric subunit is GS 

/ 215 -104.1 5.1 \ 

HRenger = -104.1 220 32.6 . (15) 

\ 5.1 32.6 0 J 

The eigenenergies of the system ar given by = 

322.85 cm~^,E_ = 119.13 E^ = —6.98 cm~^ 

which yields the eigenperiods = (27r/i/A£’y) : 

T_|_= 0.163ps, T _3 = O.lOOps, T +3 = 0.264ps. Due 

its structure, the trimer is a chain composed by a pair of 
chromophores ( 1 , 2 ), degenerate in energy and forming a 
strongly coupled dimer, and a third chromophore moderately 
coupled with the second one only. Since in the following 
we suppose that the exciton starts from site 1 , we expect a 
prominent role of the dimer in the dynamics, at least in the 
first tens of femtoseconds. 


g{t) = E^\g{0)). (13) 

In the notation above, the propagator has the form ICt't = 
gi£(t -t) q-jjg efficiency of the transport can be evaluated as 
the leak of the population Pe{t) = (e|p|e) of the exit site e 
towards the sink: 

t 

ri{t) = 2ktrap J {G\p\e). (14) 

0 

The overall efficiency of the process is obtained by letting t 
oo. 

While its Markovianity limits the faithful description of de¬ 
coherence processes actually taking place in real photosyn¬ 
thetic systems, the model retains the basic and commonly ac¬ 
cepted aspects of decoherence, that acts in the site basis: al¬ 
beit in a complex non-Markovian way, the protein enviroment 
measures the system locally (i.e., on each site), thus destroy¬ 
ing the coherence in the site basis and creating it in the ex¬ 
citon basis. Note that the formalism can also be applied to a 
‘dressed’ or polaronic basis where we include strong interac¬ 
tions between chromophores and vibrational modes. That is, 
to apply the DH method, one only needs a model in which an 
exciton hops between sites, dressed or undressed. The model 
is therefore suitable to readily implement the decoherent his¬ 
tories paradigm and to spot the main basic features we are in¬ 
terested in and that are at the basis of the success of ENAQT. 

The FMO unit has 7 chromophores and a complex energy 
and coupling landscape with no symmetries. Energies and 
couplings (i.e., the Hamiltonian H) can be obtained by differ¬ 
ent techniques: They can be extracted by means of 2D spec¬ 
troscopy as in El or computed through ab initio calculations 
as in 12511 ■ with similar but not exactly equal results. This very 
complex struxture makes EMO far from ideal as a first exam¬ 
ple to study. We thus prefer to start by working with a much 
simpler, yet fully relevant subsystem: the trimeric unit com¬ 
posed by the sites 1, 2 and 3 of the EMO complex in the nota¬ 
tion of 11^ 1^ ). The first chromophore is the site in which the 


In order to show how the DH analysis can be implemented, 
in the following we are going to consider histories in the site 
and the energy bases, with N projections at times nAt, n = 
1, ..,N. We first use the coherence function C{P, N, At) in¬ 
troduced above ( fTOb to evaluate the global coherence of the 
exciton transport process. In order to test the behavior of C 
for different values of dephasing, in Eig[T]we first plot C as a 
function of the time interval At between projections for two 
values of the dephasing rate: t) 7 = 0 , corresponding to the 
full quantum regime (EiglU) for the site basis (a); ii) 7 = 10 
corresponding to an intermediate value of dephasing (Eig. [T]l 
for the site basis (b) and the energy basis (c). 

Before entering the discussion of the various regimes, we 
note that as a function of the number of projections N all 
curves display the expected behavior: The increase (decrease) 
of the number of projections corresponds to a temporal fine- 
graining (coarse-graining) of the evolution; therefore, an in¬ 
crease (decrease) of N should imply an increase (decrease) 
of the amount of coherence between histories. As shown in 
Eig. ([TJ the function C correctly reproduces the fine (coarse) 
graining feature: the qualitative behavior of C as a function 
of At is not affected by the choice of N, while an increase 
of N corresponds, at fixed At, to an increase of C. We will 
therefore use in the following the value N = A that allows 
for a neat description of the phenomena and for a reasonable 
computational time. 

As for the behavior at fixed N, we have that in the full quan¬ 
tum regime (7 = 0 ), the system obviously displays coherence 
in the site basis only since 

Tx[\E,){E,\e-^^^^pE^^^\E,){E,W = {E,\p\E,)5,^, 

and the decoherence matrix V in the energy basis is diago¬ 
nal and independent on At and N. This simply means that 
in the full quantum regime histories in the exciton basis are 
fully decohered, since the system is not able to create coher¬ 
ence among excitons. Still in the full quantum regime, in the 
site basis, the coherence oscillates as the exciton, starting at 


7 




(a) 


(b) 



C10 



A t (ps) 


(c) 


(d) 


Figure 1: Coherence function C[P, N, At) as a function of At for the trimer with HamiltonianllSIevaluated: (a) in the site basis for 7 = 0 , 
(b) in the site basis for 7 = 10ps~^ (c) in the exciton basis for 7 = 10ps“^ (d) in the exciton basis for different values of dephasing 7 . 


site 1 , goes back an forth along the trimer, and the evolution 
builds up coherence in this basis, see FiglT](a). In this regime, 
the trimer can be approximately seen as a dimer composed by 
the hrst two chromophores, and the exciton performs Rabi os¬ 
cillations with a period given by T_|_= 160 /s; C oscillates 

with half the period; for At = 80 fs the exciton is migrated 
mostly on site 2 and C has a minumum - which is different 
from zero since the exciton is partly delocalized on site 3, and 
the system therefore exhibits a non vanishing coherence. 

For intermediate values of 7 ss 10ps“^ , Fig[T](b), the co¬ 
herence in site basis as measured by C correctly drops down 
at At > 7 “^|| 3 ]. The dephasing has a strong and obvious ef¬ 
fect on the coherence between pathways: coherence in this 
basis is a monotonically decreasing function of 7 . This is 
well highlighted by the global coherence function C, whose 
maximal values are reduced by a factor of ^ 3 with respect 
to those corresponding to full quantum regime. After a time 
Tdecoh = 7 ”^ the histories are fully decohered. Indeed, due 
to the specific model of decoherence (fTTT) . which amounts to 
projective measurements on |i)(f| at each site with a rate 7 , 
the system kills the coherence in the site basis, which in turn 
corresponds to the stable pointer basis for this model 101 , i.e., 
the basis in which the density matrix is forced to be diagonal 
by the specihc decoherence model. On the other hand, and for 
the same reason, the dynamics starts to build up coherence in 


the exciton basis \Ei){Ei\, see Fig. [Uc). However, this co¬ 
herence is later destroyed - on a time scale of approximately 
0.2 ps - since the stationary state of the model is the identity. 
This effect is even more evident if one compares the behavior 
of C in the exciton basis for different values of 7 , as shown in 
Fig. Ed): C grows with 7 and it lasts over longer time scales. 
This feature is coherent with the expectations: the equilib¬ 
rium state for high 7 is the identity. Due to the projections 
implemented by the environment in the site basis, the system 
is forced to create coherence in the exciton basis. When 7 
is very high a quantum Zeno effect takes in, the dynamics is 
blocked, and the time required to reach the equilibrium, and 
to destroy coherences in all bases, consequently grows. 

This first analysis therefore shows that C is indeed a good 
candidate for assessing the global coherence properties of 
quantum evolutions. For a fixed number of projections N, 
C{At) can be interpreted as a measure of the global coher¬ 
ence exhibited by the dynamics over the time scale At. 

We now analyze in detail the specihc features of quan¬ 
tum transport for the trimer. The dynamics starts at site 1 
and evolves by delocalizing the exciton on the other chro¬ 
mophores. In order to study this process, we hrst use a mea¬ 
sure of delocalization introduced in j^ for the study of LHCII 
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complex dynamics: 

Hit) = - Y^Piit) Inpj(t) ( 16 ) 

i 

that is simply the Shannon entropy of pi{t), the populations 
of the three chromophores. This measure allows one to follow 
how much the exciton gets delocalized over the trimer with 
time and in different dephasing situations: 0 < (f) < ln(3), 

1. e., TL is zero when the exciton is localized on a chromophore 
and it takes its maximal value ln(3) when the population of 
the three sites are equal. In Fig. |2]we plot both 'H(f) and 
the population P3{t) of site 3 for different values of 7 . Due 
to the presence of interference, in the mainly quantum regime 
(7 = 0 .1,1 ps~^), the exiton first delocalizes mainly over the 
dimer and partly on the third site: The first maximum corre¬ 
sponds to f = 40 /s = 1/4 T+- when the system builds up a 
(close to uniform) coherent superposition between sites 1 and 

2 , while a non negligible part of the exciton is found in site 

3; indeed = 40 fs) Ri 0.75 > 0.69, the last value corre¬ 
sponding to ln(2) i.e., to a uniform superposition over the sites 
1 and 2 only. As the dynamics of the sytems extends to later 
times we see that 77 (f) and p^ {t) have an oscillatory behavior, 
whose main period is 1/2 T_|_, and which approximately cor¬ 

responds to Rabi oscillations between site 1 and 2, although 
the initial state fully localized in site 1 cannot be rebuilt due to 
the presence of site 3. As for the transport, we see that in this 
regime the system cannot take advantage of the initial fast and 
high delocalization: the exciton bounces back and forth over 
the trimer. In the intermediate regime 7 = 16 ps~^, due, as 
we will later see, to the selective suppression of interference 
processes, the initial speed up in delocalization is sustained by 
the dynamical evolution, and the transfer rate to site 3 is cor¬ 
respondingly increased. For very high values of decoherence 
(7 = 100) the role of initial interference is suppressed and 
the initial speed-up disappears: the environment measures the 
system in site basis at high rates and the delocalization process 
is highly reduced. 

The optimal delocalization occurs in correspondence of 
7 Ri 16 and it can be interpolated with a double ex¬ 
ponential function 


%-y=lG(t) = Co + CiC + C2e 

with Co = 1.098 = In 3, a = -0.84, C 2 = -0.373. The 
first time scale ti = 23 fs describes the initial fast quantum 
delocalization process described above, while the second time 
scale T 2 — 238 fs the slower subsequent delocalization and 
the reaching of the equilibrium situation, ?{(t = 00) = ln(3). 

We now pass to systematically analyze the behavior of 
the coherence of the evolution with respect to the strength 
of the interaction with the environment and its relevance for 
the energy transport process. As a first step we plot both 
77 (r = NAt) and C{At) for different values of 7 , Fig. Q. 
The plots show that the coherence function exhibits the re¬ 
quired behavior: For small 7 = 0.1 , C{At) oscillates with 
period 1/2 T_|_, following the Rabi oscillations of the dimer. 



(a) 



(b) 

Figure 2: For the trimer (E): (a) Delocalization 77(i) and (b) P3{t) 
population of site 3 as a function of time for different values of 7 = 
0.1,1,16,100 ps-i 



C,y=1 ps 
H,y=1 ps 
C, y=16 ps 
H, y=16 ps 
C, y=60 ps 
H, y=60 ps' 


1 

1 

1 

1 

1 

1 


Figure 3: Delocalization 77 (r = NAt) and coherence function 
C{At) for 7 = 1 ,16, 60 ps~^ 


The minima occur at n/4 T_|_, showing that the exciton is 

“partially” localized on site 1 or 2 , and partially delocalized 
on site 3. As 7 grows, the system becomes unable to create 
coherence on large time scales; the decay of C{At) is mir¬ 
rored by a the reduction of the amplitude in the oscillations of 
77(r = NAt). 

We now focus on the relevant time scales for the initial 
fast delocalization process highlighted by our previous anal- 
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Figure 4: Average coherence of the evolution Qt^ ( 7 ) for the trimer 
G3 for different values of Td 


ysis, which are of the order of tens to hundreds of femtosec¬ 
onds. We therefore introduce the following average measure 
of global coherence of the evolution 

QtAi) = —[ C{At)dAt. 

Td Jo 

Q-Tdil) is the average of the coherence exhibited by the dy¬ 
namics of the system at the time scales At G (0,rd). In 
Fig. |4] we show Qrd ( 7 ) for the trimer (flST l in the site basis 
for different values of Td- We first focus on the behavior of 
Qrdil) for values of dephasing in the range 7 G ( 0 , 1 ) ps~^. 
In this range, for small time scales Td = 20 to 200 fs the 
average global coherence Qr^il) is approximately constant 
and equals the value attained in the full quantum regime i.e., 
QrAl) ~ QtAi = 0 . 1 ) . For larger time scales (r^ Ri Ips) 
( 7 ) rapidly decreases with 7 . This analysis shows that the 
behavior of Qr^ ( 7 ) matches the expectations: the higher 7 the 
smaller the time scales over which decoherence takes place, 
the lower the global coherence of the dynamics. Along with 
C{At) the functional Qr^il) is therefore in general a good 
candidate for the evaluation of the global coherence of open 
quantum systems evolution. As for the transport dynamics, we 
focus on the timescale identihed with the analysis of H(f) for 
optimal dephasing; for Td = ti = 20 ps and Td = 40 ps we 
see that the system indeed retains most of the average coher¬ 
ence of the purely quantum regime up to the optimal values of 
decoherence (7 = 16ps“^ in the hgure), losing it afterwards; 
this is a clear indication that this phenomenon is at the basis 
of the the fast initial delocalization process. Over longer time 
scales, the relevance of coherence is highly suppressed. 

We now deepen our analysis about the relevance of the co¬ 
herence of the evolution for the energy tranfer efficiency. To 
this aim we focus on the basic feature that distingushes the 
classical and the quantum regime: interference. In particular 
we focus on the sub-block X >3 of the decoherence matrix V 
pertaining to the third chromophore, which describes the set 
of histories in site basis ending at site 3. Due to interference 
the probability of occupation of the site 3 at time r = NAt 
can be written in terms of the the histories ending at site 3 
P 3 (r) = wo{t) -f I 3 (t), see (gji . 



7 = 0.1 

7=1 

7=10 

7=100 


ps 


ps 


ps 



Y=5 ps 
Y=8 ps 
Y=11 ps 
Y=14 ps 
Y=18 ps 
Y=30 ps' 


1 

1 

1 

1 

1 

1 


(b) 


Figure 5: (a) interference ^^{t) of histories ending in site 3 for the 
trimer as a function of r = NAt for different values of 7 = 
0 .1, 1,10,100; (b) ^^{t) for different intermediate values of 7 


In Fig. |5]we show Is^t) for different values of dephasing. 
One has different regimes: for 7 ^ 1 , the set of histories in 
site basis is fully decohered; I 3 (t) ~ 0 , the histories do not 
interfere with each other and P 3 (r) r; W 3 {t) i.e, the proba¬ 
bility is simply the sum of diagonal elements of D 3 . In the 
mainly quantum regime 7 < lps“^, P 3 (t) f W 3 {t): af¬ 
ter the initial positive peak the histories interfere with each 
other, globally the interference is mostly negative and there¬ 
fore P 3 {t) < W 3 {t). For intermediate values of decoherence 
7 “^ R 10 ps the interference has a positive peak and then 
reduces to zero. While the hrst initial hngersnap of positive 
interference that takes place in the hrst r 80 fs is common 
for all curves corresponding to small and intermediate values 
of 7 , the main effect of the bath is displayed after this initial 
period of time: the decoherence gradually suppresses inter¬ 
ference, both the positive and the negative one', however, for 
intermediate values of y the effect is stronger as for the neg¬ 
ative part of the interference patterns. The environment thus 
implements what can be called a quantum recoil avoiding ef¬ 
fect: it prevents the part of the exciton that - thanks to con¬ 
structive interference - has delocalized on site 3 to How back 
to the the other sites. 

In order to evaluate a possible advantage provided by the 
initial speed-up in the delocalization process and by the inter- 
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Figure 6: For the trimer 01 (a) Average positive (green), negative 
(blue), total (red) interference (Is) of histories ending in site 3 as a 
function of 7 ; (b) total average interference for site i — I 

(red), i = 2 (blue), i = 3 (green) 7 . 


ference phenomena showed above one has to take into account 
another relevant time scale of the transport process; the trap¬ 
ping time. Indeed, if the system is to take advantage of the 
fast delocalization due to the coherent behavior, the exit of 
the exciton should take place on time scales of the order of 
the delocalization process. The theoretical and experimen¬ 
tal evidences show that this is the case; the trapping time 
Ttrap — '^^6 FMO complex is estimated in the lit¬ 

erature to be of the order of 0.2 ps i.e., the exit of the exiton 
starts soon after the fast delocalization due to quantum coher¬ 
ence has taken place. The role of the interference between 
paths, in particular those leading to site 3, can therefore be 
appreciated by numerically evaluating 

(If) = / If (r)dr (17) 

'trap Jo 

i.e., the average over the trapping time scale of Ttrap = 200/s 
of the total (/3 = T ot), negative (/3 = —) and positive (J3 = -F) 
average interference between the histories ending in site i, 
with = (if") -F (Ij~)- In particular, in Fig. IS a) the 

different kinds of interference are plotted for histories termi¬ 
nating at site 3; on average, the negative interference highly 
reduces the total interference for small values of decoherence 
strength; when 7 Ri 10ps“^, (I^) vanishes, the average total 
interference equals the positive one (7^°*) = (if^), and it is 
maximal for values of 7 comparable to those that maximize 


ri, y =0 ps'^ - 

W 3 , y=0 ps'^ . 

I3, 1=0 ps'^ . 

U, y=16 ps'^ 

W 3 , y=16 ps'^ 

I 3 , y=16ps'^ 

U, y=60 ps'^ - 

W 3 , y=60 ps'^ . 

I 3 , y=60 ps'^ . 


0 1 2 3 4 5 

t (ps) 

Figure 7; For the trimer d: transport efficiency 17 (f), integrated 
weight Wiit) = 2ktrap Jq W 2 ,{T)dr and integrated interference 
73 (f) = 2ktrap Jq lair)dr for pathways ending at site 3 for dif¬ 
ferent values of 7 = 0 . 1 , 16, 60 ps~^ and ktrap = 5ps“^ 



77(f) (« 16 ps ^). In Fig. | 6 l(b), we compare the behavior of 
(l7°‘) for all sites. The results again suggest that decoher¬ 
ence acts on the interference provided by the quantum engine 
in order to favor the flow of the exciton towards the exit chro- 
mophore; the average positive interference between histories 
ending at sites 2 and 3 grows in modulus with 7 and attains a 
maximum for intermediate values of decoherence; while the 
average negative interference between histories ending at site 
1 decreases and attains a minimum for intermediate values of 
7 . The combined effect of decoherence and interference thus 
helps depopulating site 1 and populating site 2 and 3. 

We can now tackle one of the most relevant aspects of our 
discussion; the net effect of the above described phenomena 
on the overall efficiency of the transport. The latter can be 
fully appreciated by evaluating the efficiency of the process 
(fT4l) and by recognizing that, in the decoherent histories lan¬ 
guage, it can be expressed as; 


t 

ViO = ‘Iktrap J P 3 {T)dT = Wsit) + 73 ( 7 ) 

0 

where t = NAt and Ws^t) = 2ktrap Jq W3{T)dT, Is^t) = 
2ktrap Jq l 3 {T)dT. This split allows one to appriciate the role 
of interference for the efficiency. In Fig. (|7]) 77 is plotted for 
different values of dephasing. In agreement with what dis¬ 
cussed above, we have three regimes; for very small values of 
7 the overall efficiency is poor; this is due to the presence of 
high negative interference that in average prevents the exciton 
to migrate to the exit site. For large values of 7 the inter¬ 
ference processes are completely washed out and the system 
cannot take advantage of the fast quantum delocalization. For 
intermediate (optimal) values of 7 only the negative interfer¬ 
ence has been washed out; l 3 (r) is positive, it acts on short 
time scales, and it provides on average an enhancement of the 
global efficiency. 
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These results, within the limits of the simple model of de¬ 
coherence taken into account, undoubtedly show for the first 
time that the so called ENAQT phenomenon can well and 
properly be understood both qualitatively and quantitatively 
within the decoherence histories approach, i.e., in terms of 
very the basic concepts of coherence and interference between 
histories. The often recalled “convergence” of time scales 
or “Goldilocks” effect ( 1^ 1 in biological quantum transport 
systems seems therefore to be well rooted in the processes dis¬ 
cussed above; if decoherence is too small the system shows 
both positive and negative interference (see Fig. Ell, the delo¬ 
calization has an ocillatory behavior, and the exciton bounces 
back and forth along the network thus preventing its efficient 
extraction. If instead decoherence is very high one has that 
the complete washing out of intereference and coherence im¬ 
plies the delocalization process to be very slow, no matter how 
fast the trapping mechanism try to suck the exciton out of the 
system. In order to take advantage of the effects of quantum 
coherent dynamics: i) the bath must act on the typical time 
scales of quantum evolution in order to implement the quan¬ 
tum recoil avoiding process; ii) the extraction of the exciton 
from the complex, characterized by ktrap, must then start soon 
after the initial fast delocalization has taken place. Should 
the extraction take place on longer time scales, the benefits 
of the fast initial delocalization would be spoiled; Waiting 
long enough, the system would eventually reach together with 
equilibrium a decent delocalization even for moderately high 
values of 7, but in this case the transfer would be obviously 
much slower. 


V. FMO 

The above arguments can be easily applied to the whole 
FMO complex. Fig. [ 8 ] and |9] show the application of the de- 
coherent histories method to excitonic transport in FMO. The 
main features of the behavior of C, Q,l 3 and 77 are main¬ 
tained although obvious differences can be found since the dy¬ 
namics in now determined by the interplay of different eigen- 
periods and interference paths are more complex. In particu¬ 
lar, Fig| 8 jb), one can observe a revival of positive interference 
I 3 for small values of 7 ~ 1 , that does enhance the efficiency 
for t Ri 1.5ps, Fig. lUa); but this is not sufficient to com¬ 
pensate the initial and subsequent negative interference, thus 
impeding the reach of optimal values of 77 . In general, com¬ 
pared to the trimer and as suggested by Fig. |9ja) the maximum 
average positive coherence on short time scales is attained for 
smaller values of 7. The overall picture is not significantly af¬ 
fected if one decides to start the dynamics from site 6 instead 
of site 1 , as it often is reported in the literature. 


VI. CONCLUSIONS 

The decoherent histories approach provides a general the¬ 
ory to study the distinctive feature exhibited by quantum sys¬ 



la) 



(b) 



(c) 


Figure 8: FMO: (a) C{At) for different values of 7; (b) l 3 (At) 
for different values of 7; (c) 'H{t = NAt),C{At) for 7 = 
1,16, 60 


terns: coherence. However, despite its generality and foun¬ 
dational character, in order to measure the effects of coher¬ 
ence and decoherence the DH approach needs to be comple¬ 
mented with a quantitative way to condense the information 
contained in the basic object of the theory, i.e., the decoher¬ 
ence matrix V. In this paper we introduce a set of tools that 
allow one to assess the (global) coherence properties of quan¬ 
tum (Markovian) evolution and that can be used to relate the 
coherence content of a general quantum dynamical process to 
the relevant figure of merits of the given problem. We first 
define the coherence functional C{P, N, At), that can be in- 
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(a) 



W, y=0. 

I, y=0 . 

Ti, y=16 
W, y=16 
1 , 1^16 

Ti, y=60 - 

W, 1^60. 

I, y=60 . 


Figure 9 : FMO: (a) average interference < Is > for ktrap = 
5 ps~^\ (b) transport efficiency integrated weight W3{t) and 
integrated interference /3(f) for pathways ending at site 3 for differ¬ 
ent values of 7 = 0,1, 16 ps“^ and ktrap = 5 ps~^ 


terpreted as a measure of the global coherence exhibited by 
the dynamics in the basis P over the time scale At. While 
this measure is completely general, one can further introduce 
other relevant tools tailored to the specific system and type 
of system-environment interaction at hand. We thus focus on 
a simple yet paradigmatic model of environmentally assisted 
energy transfer where coherence effects have been shown to 
play a significant role in determining the efficiency of the pro¬ 
cess: a trimeric subunit of the Fenna-Matthews-Olson photo¬ 
synthetic complex. Based on 1) and C{P, N, At) we define: 
a) a measure Qri'j) able to characterize the average coher¬ 
ence exhibited by the dynamics of the system over the time 
scales At G (0, t ) for a fixed value of the dephasing 7 ; b) 
a measure of the average interference < It > occuri'ing be¬ 
tween the histories ending at a given “site” i. 

Within the specific model, we first thoroughly assess the 
consistency of the behavior of C{P,N,At) in the various 
regimes. We then show how the introduced tools allow to 
study the intricate connections between the efficiency of the 
transport process and the coherence properties of the dynam¬ 
ics. In particular we show that the delocalization of the exci- 
ton over the chromophoric subunit is strongly affected by the 
amount of (average) coherence allowed by the interaction with 


the bath in the first tens to hundreds of femtoseconds. If the 
system-bath interaction is too strong, coherence is suppressed 
alongside the interference between different histories, in par¬ 
ticular those ending at the site where the excitation leaves the 
complex. If the interaction is too weak the system exhibits 
high values of coherence even on long time scales, but it also 
exhibits negative interference between pathways ending at the 
exit site, a manifestation of the fact that the exciton bounces 
back and forth over the network thus preventing its efficient 
extraction. In the intermediate regime i.e., when the different 
time scales of the system (quantum oscillations, decoherence 
and trapping rate) converge, the system shows high values of 
coherence on those time scales. The action of the bath has a 
quantum recoil avoiding effect on the dynamics of the exci¬ 
tation: the benefits of the fast initial quantum delocalization 
of the exciton over the network are preserved and sustained 
in time by the dynamics', in terms of pathways leading to the 
exit site, the action is to selectively kill the negative interfer¬ 
ence between pathways, while retaining the initial positive 
one. These effects can be explicitly connected to the over¬ 
all efficiency of the environment-assisted quantum transport: 
the gain in efficiency for intermediate (optimal) values of de¬ 
coherence can thus be traced back to the basic concepts of 
coherence and interference between pathways as expressed in 
the decoherent histories language. 

While the specific decoherence model used (Haken-Strobl) 
is an oversimplified description of the actual dynamics tak¬ 
ing place in real systems, we believe that our analysis al¬ 
lows to spot out the essential features that may determine the 
high efficiency of the transport even in more complex system- 
environment scenarios. 

In conclusion, the tools introduced in this paper allow to thor¬ 
oughly assess the coherence properties of quantum evolutions 
and can be applied to a large variety of quantum systems, the 
only limits being the restriction to Markovian dynamics and 
the computational efforts required for high dimensional sys¬ 
tems. However, the extension to non-Markovian realms is 
indeed possible lO . and the use of parallel computing may 
allow the treatement of reasonably large systems. 
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